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■ Abstract: We investigate the production of a heavy charged Higgs boson at hadron 

colliders within the context of the MSSM. A detailed study is performed for all important 
production modes and basic background processes for the ttbb signature. In our analysis 
we include effects of initial and final state showering, hadronization, and principal detector 
effects. For the signal production rate we include the leading SUSY quantum effects at 
high tan/3 > mt/mt,. Based on the obtained efficiencies for the signal and background 
we estimate the discovery and exclusion mass limits of the charged Higgs boson at high 
values of tan (3. At the upgraded Tevatron the discovery of a heavy charged Higgs boson 
(Mjj± > 200 GeV) is impossible for the tree-level cross-section values. However, if QCD 
and SUSY effects happen to reinforce mutually, there are indeed regions of the MSSM 
parameter space which could provide 3 a evidence and, at best, 5 a charged Higgs boson 
discovery at the Tevatron for masses M H ± < 300 GeV and M H ± < 250 GeV, respectively, 
even assuming squark and gluino masses in the (500 — 1000) GeV range. On the other 
hand, at the LHC one can discover a as heavy as 1 TeV at the canonical confidence 
level of 5 a; or else exclude its existence at 95% C.L. up to masses ~ 1.5 TeV. Again the 
presence of SUSY quantum effects can be very important here as they may shift the LHC 
limits by a few hundred GeV. 
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1. Introduction 



The full experimental confirmation of the Standard Model (SM) is still waiting for the 
finding of the Higgs boson. Last LEP results, suggesting a light Higgs boson of about 
115 GeV [Q], are encouraging, but we will have to wait new data from the upgraded Tevatron 
or the advent of the LHC to see this result either confirmed or dismissed g, ||, (|, [Bj. 
Moreover, at the LHC a large amount of top-quark pairs will be produced, allowing for 
high-precision measurements of all its properties, providing strong checks of the SM, or 
new physics signals Q. Even if a neutral Higgs boson is discovered, the principal question 
will still be present at the forefront of Elementary Particle research, namely: is the SM 
realized in nature or does a model beyond the SM take place with extended Higgs sectors 
of various kinds (extra doublets, singlets, even triplets)? In most of these extensions, the 
physical spectrum contains charged Higgs bosons and this introduces a distinctive feature. 
For example, in the general two-Higgs-doublet model (2HDM) Q one just adds up another 
doublet of scalars and then the spectrum of the model contains three neutral Higgs bosons 
and two charged ones, the latter being commonly denoted by H^. This is also the case of 
the Higgs sector of the Minimal Supersymmetric Standard Model (MSSM) [||, which is a 
prominent Type II 2HDM of a very restricted kind. 

While the detection of a charged Higgs boson would still leave a lot of questions 
unanswered, it would immediately offer (in contrast to the detection of a neutral one) 
indisputable evidence of physics beyond the SM. Then the next step in charged Higgs boson 
physics would be the precise measurement of its mass and couplings to fermions, to give us 
the understanding of whether those couplings and masses are compatible with the ones in 
the MSSM or belong to a more general 2HDM. Of course this task cannot be accomplished 
without including the information provided by radiative corrections. These are not only 
potentially large in the computation of the MSSM Higgs boson masses themselves Q but 
also in the interaction vertices and self-energies. The latter features have been studied in 
great detail for the decay of the top quark into a charged Higgs boson [|i]] and also for the 
hadronic decays of the Higgs bosons of the MSSM [|ll], 12]. 



Focusing on the charged Higgs boson, the information we have nowadays is rather 
limited. First, there is a direct LEP limit, M H ± > 79.9 GeV, for a 2HDM #± decaying 
exclusively into tv t |l~3[j . Second, Tevatron direct and indirect searches for a place 
constraints on the M H +-BR{t — > bH + ) plane [14, [0|], which are usually translated to 



the M H + -tan/3 plane once the relevant MSSM parameters are fixed p5| , 16, 17] ; here 
tan/3 stands for the ratio between the two Higgs doublets vacuum expectation values, 
tan/3 = Vz/vi @. Such searches have been done for M H + < mt — rrib so far. Finally, 
the bound coming from the rare process b — * S7, for a pure Type II 2HDM, places a tight 
constraint in the M H + -tan/3 plane, although in the MSSM lighter masses can be traded 



for a constraint in the A±—\i plane when the -R-parity is conserved |18| or broken [19]. This 



constraint has been found to be robust under higher order effects [20]. Unlike the direct 
searches at LEP and the Tevatron, which rely on the electromagnetic coupling and on 
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a few, quite universal, decay modes, the weak point of the indirect limits is their model 
dependence. 

Several other processes have been proposed in the literature that would allow detection 
of a charged Higgs boson in various kinematic ranges. Among others, there are: pair 
production at e + e~ [^] and hadron [22] machines; associated production with a W boson 



in e + e p3] and hadron [24, 25] colliders. Indirect searches in B-physics observables 



(apart from b — > S7) are also promising, specially in b — > sl + l~ and b — > ctu t [26]. 
Very high precision measurements on r and K physics can give also interesting indirect 
limits p7| ]. Further production channels in hadron colliders specific to SUSY models include 
the possibility of the charged Higgs boson being produced in cascade decays of strong 
interacting SUSY particles, as well as associated production with squarks |28|]. These 
channels are mostly useful in intermediate ranges of tan/?, and are therefore complementary 
to the channel studied in the present work. 

The aim of this paper is to extend the analysis of the associated production with 
top (and bottom) quarks at hadron colliders pp(pp) — ► H + tb + X to both the Tevatron 
and the LHC. In this work we will be focusing exclusively on the MSSM case. The general 
2HDM case will be addressed elsewhere ]^9| ; following [^0| we also expect here sizeable 
quantum effects, though of very different origin and distinguishable from the MSSM ones. 

From our point of view, the relevance of the process pp(pp) — > H + ib + X is manifold. 
A positive signal would be instant evidence of new physics and at the same time it would 
strongly indicate a large value of tan (3. In fact, as we shall see, the great virtue of this pro- 
cess is its ability to test the charged Higgs boson coupling to the third generation of quarks 
and leptons, which in general 2HDM's (and in particular in the MSSM) can be greatly 
augmented (or suppressed) not only at the tree-level but also due to model-dependent 
quantum effects. Such an enhancement could be critical for heavy production, espe- 
cially for the Tevatron. Furthermore, if one would be able to correlate the quantum effects 
on the charged Higgs boson Yukawa coupling with the corresponding effect in the measured 
value of the neutral Higgs boson Yukawa couplings, it would mean a strong hint for a MSSM 
Higgs sector, even if no supersymmetric particle is detected at all [11, ll^. |3l], |32fl. That this 



kind of scenario is possible is corroborated by the fact that the process pp(pp) — ► H + ib + X 
is affected by exceptional SUSY quantum effects, namely effects that do not necessarily 
vanish in the limit of large MgusY masses [^] . This remarkable property has been greatly 
emphasized from the phenomenological point of view for applications in both high- and 
low-energy processes in [1C, [H|, 12,31, |32], 34, and |36|] respectively. 1 Certainly, here is 



another place where it might play a fundamental role to help uncovering the SUSY nature 
of the charged Higgs boson in high-energy experiments. Indeed, we will show that these 
Higgs boson production processes could receive large radiative corrections even if the di- 
rect cascade SUSY channels mentioned above are kinematically closed. In other words, the 
Higgs production modes addressed in the present study are not only generally competitive 
with the direct SUSY channels, but they could also be highly relevant by themselves. 



For a detailed analysis of the neutral MSSM Higgs bosons production at the Tevatron, with the inclusion 
of leading SUSY effects, see @. 
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At this point, a review of the previous works on pp(pp) — > H + ib + X is in order. To 
the best of our knowledge, the first analyses on this process in the literature are the works 
of Refs. [3J|. Recently, there have been new contributions retaking this issue in more 
detail |4(], ^j], [42], |43|, ^4|] . Nevertheless all of these works stick to a tree-level computation, 
in spite of the fact that some of them explicitly admit that the sort of charged Higgs boson 
they are dealing with is of the MSSM type. Therefore, they are unavoidably affected by 
some drawbacks. On the one hand, the tree-level approach can be quite inaccurate (e.g. 
gluino radiative corrections are potentially of order 1 for large tan/3 and have, a priori, no 
definite sign) and, on the other, it is unable to distinguish between a generic and one 
with the particular couplings of the MSSM. 

In Ref . [fij| , a first treatment of the SUSY radiative corrections was presented although 
it was not sufficiently complete and moreover no background analysis was attempted. 2 
Therefore, we believe that considerably more discussion is needed before any definite pre- 
dictions on the discovery limits for the charged Higgs boson as a function of the MSSM 
parameter space can be made. 

The paper we are presenting here gives a significant further step in this direction. 
It adds, with respect to previous studies, a simultaneous treatment of the leading SUSY 
radiative corrections (both strong and electroweak) including an analysis of the off-shell 
effects. Also, in contrast with previous works we perform a beyond-the-parton-level simu- 
lation of events, which includes the toy-detector simulation, jet fragmentation, initial and 
final radiation effects. We show that some of those effects, such as final radiation, do play a 
crucial role in the charged Higgs boson studies. Another new aspect of the present paper is 
the proper kinematical analysis of gg — > H + ib and gb —* H + i processes after their correct 
combination. 3 First results of the present analysis were presented in ||46f| . 

Last but not least, it will also be important to discuss the possible effect from the 
conventional QCD corrections, which unfortunately are not available for pp(pp) — > H + ib + 
X at present. However, two independent calculations for the related process in the Standard 
Model, pp(pp) — ► Hti+X, have recently appeared, namely the associated SM Higgs boson 
production at the Tevatron and the LHC, which have been carried out at the next-to-leading 
order (NLO) in QCD 0, |f]. The result is that at the Tevatron the NLO effects for SM 
Higgs boson production off top quarks are negative and can be approximately described by 
a if-factor ranging between 0.8— 1 whereas at the LHC they are positive and the iT-factor 
lies between 1.2 — 1.4 depending on the renormalization and factorization scales. Obviously 
this process bears relation to ours, pp(pp) — > H + tb + X, and so a few comments on the 
expected NLO QCD corrections on the latter are in order (see Section ^). 

The rest of the paper is organized as follows: in Section || we discuss the general 
procedure for the computation of the cross-section. In Section || we perform a signal and 
background study based on the PYTHIA6.1 simulations, and work out kinematical cuts 
and a suitable strategy to suppress the background and extract the signal in the most 
efficient way. Then we present signal and background efficiencies and signal rates that 

2 A close comparison of the results from Refs. (fB[ can be found in [Q. 

3 Ref. ^lj makes also the proper combination of the gg — > H + ib and gb — > H + i channels for the total 
effective cross-section, however no analysis of the differential cross-section is performed. 
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could be viable at the Tevatron and LHC. In Section || we discuss in detail the role of the 
SUSY corrections both analytically and numerically, and then evaluate their impact on the 
discovery limit or exclusion region of the charged Higgs boson. The last section, Section || 
is devoted to the summary and conclusions. 

2. Cross-section computation 

The relevant charged Higgs boson production processes under study are the following: 



pp{pp) -> H + tb + X Tevatron (LHC) (2.1) 



At the parton level, the reaction ( |2.lD proceeds through three channels: i) qq- annihilation 
for light quarks 4 



qq 



H+tb, (2.2) 



where q = u,d (the s contribution can be safely neglected), a channel only relevant for the 
Tevatron |4l], ^5|; ii) gg-fusion 

gg^H + tb, (2.3) 

which is dominant at the LHC, but it can also be important at the Tevatron for increasing 



H+ masses [|4l] , fJf|; and finally there is the iii) bottom-gluon 2-body channel 

bg -» H+i . (2.4) 

We remark that the latter can also be significant because the bottom quark mass, 
m&, is small with respect to the energy of the process, and therefore parton distribution 
functions (PDFs) for b and fc-quarks (i.e. 6-densities) have to be introduced, allowing for 
the resummation of collinear logarithms [^] . This provides an extra channel contributing 
to the cross-section which must be appropriately combined with the gg-fusion channel as 
we will comment further below. 

We will compute the cross-section for the charged Higgs boson production process (pj] ) 
at the leading order (LO) in QCD, namely at 0{a 2 s ). Feynman diagrams for the partonic 
subprocess (2.2), (|2.3| ) and (2.4) are presented at the tree-level in Fig. [l|. However, the 



QCD corrections at the next-to-leading order (NLO) or C(a|) could be important. We 
have already mentioned in the introduction that the full set of QCD corrections to the 
process of SM Higgs boson radiation off top quarks in hadron colliders, 

pp{pp) -> H ti + X , (2.5) 



have been computed at the NLO by two independent groups [47, 48]. The result is that, at 
the Tevatron, the NLO effects from QCD are negative, hence diminishing the signal cross- 
section, whereas at the LHC they are positive and so enhancing the signal. To be more 

4 We shall omit the charge-conjugate process, pp(pp) — > H~tb + X, for the sake of brevity. Including this 
process just amounts to multiplying our cross-section by a factor of 2. 
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Figure 1: Diagrams for the tree-level processes (2.3) - a), (|2.2[) - b) and (2.4) - c). 



specific, if one defines the scale no = mt+Mn /2, then the NLO effects can be approximately 
encoded in a K- factor as follows [47, ^8|. Assuming equal factorization and renormalization 
scales Q = = fj,, the i^-factor for the Tevatron ranges between K = 0.8 (for the central 
scale value \x = Hq) and K = 1 (for the threshold energy value fj, = 2 no = 2m t + Mh), 
whereas at the LHC it varies between K = 1.2 (for \x = /j,q) and K = 1.4 (for [i = 2^o)- 
These results are in agreement with the expectations from the fragmentation model of 
Ref . [50|] where an approximate calculation of the NLO cross-section can be performed in 
the limit of small Higgs boson masses. In particular, the negative sign of the QCD effects 
at the Tevatron can be understood from the dominance of the qq — > Htt partonic mode at 
the Tevatron energies (or, equivalently, for heavy quark masses such as the top quark mass) 
which is subject to (moderate) negative QCD effects both near and above the threshold. 
Clearly, the value of the K-factor can be important, and even critical, especially at the 
Tevatron where the size of the signal is not too conspicuous to allow for a comfortable 
separation of it from the QCD background. Remarkably, in the MSSM framework the 
neutral Higgs bosons production processes in association with bottom quarks can be, in 
contrast to the SM case, even more important than the associated production with top 
quarks. It can proceed through 



pp(pp) — ► hbb + X 



(2.6) 



with h = h°, H°,A° any of the neutral Higgs bosons of the MSSM These processes are 
certainly related, in fact they are complementary, to our charged Higgs boson case ([2.1|), 
and different aspects of them have already been analyzed in the SUSY context in several 
places of the literature J2), |37], 51 1. As in the charged Higgs boson case, the QCD NLO 
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corrections for the processes (|2.6| ) have not been computed yet 5 . Nevertheless we wish to 
remark that in all of these MSSM Higgs boson production processes in hadron colliders, 
with at least one 6-quark in the final state, we expect important QCD effects at the NLO 
(and even at higher orders) ||52|| . In fact, in contradistinction to the SM case (|2.5|), the 
QCD corrections for these processes should entail a large and positive if-factor for both 
the Tevatron and the LHC (j53|. This circumstance could be highly favorable to enhance 
the MSSM Higgs boson production and perhaps make one of these particles visible already 
at the Tevatron. 

In the absence of a fully-fledged calculation, we can still have some hints on the approxi- 
mate value of the QCD corrections as follows. First, the result of the NLO QCD corrections 



on the SM process (|2.5D from Refs. |47|, |48|] does not directly apply to the charged Higgs 
boson case under study fl2.1|), since the presence of a (nearly) massless quark in the final 
and intermediate states will introduce additional logarithms of the light quark mass. Sec- 
ond, part of these logarithms is taken into account by the proper combination (see below) 
of the processes fl2.3| ) and Q2.4| ), which should absorb most of the collinear logarithms due 
to the light bottom quark mass into the PDF of the bottom quark. In this respect we 
warn the reader that the convenience of the use of bottom-quark PDFs to appropriately 
describe Higgs boson production is not entirely settled down yet [53, 55 1. This is specially 
so in the case of neutral Higgs boson production (|2.6|) , where the result of adding up the 
2^2 with the 2^3 process yields a cross-section an order of magnitude larger than the 
leading order 2 — > 3 alone [§, |5q] , thus giving a strong hint that the bottom-quark PDF 
description does not approximate well the full result. Third, for the charged Higgs boson 
production cross-section ( |2.1| ) the situation is different, since the result of the combination 
of processes ( p.3[) and ( |2.4[) gives a cross-section of the same order of magnitude than that 
of process (^3|) alone - see e.g. Tables [j] and ^ below. So it is fair to think that the 
combination of subprocesses ([0|) and ( |2.4|) gives a better approximation to the charged 
Higgs boson production process than the computation based on using ( [2. 3D only. Still, 
one should keep in mind the situation of the neutral Higgs bosons, whose resolution might 
lead to a new and interesting description of these processes 53 ] . Aside from the collinear 
logarithms, other contributions to the NLO corrections exist. For example, in Ref. [54] 
the standard NLO QCD corrections to the subprocess (^4) at the LHC are computed, 
obtaining a large K-factor between ~ 1.6 and ~ 1.8 for tan/3 > 20. Taking into account 
these considerations, a large K- factor for the full process (2.1) coming from the standard 
QCD corrections (e.g. i^Q CD ~ 1.5) is not ruled out. 

On account of the previous discussion, we will make use of pole quark masses through- 
out our study. The use of the running quark mass would be justified only if we would have 
good control on the value of the remaining QCD corrections. It is well known that in decay 
processes, like for instance t — > bH + , the use of the bottom quark running mass Trib(Q) 
actually accounts for most of the QCD virtual effects (at Q = mt) J56| |. But in production 
processes the QCD corrections cannot be parametrized in this way. Indeed, a clear hint 



5 Notice that in the MSSM the process (|2.E|), with H replaced with any of h = h° , H°,A°, is not favored 
at high tan and moreover for tan /3 — 0(1) the QCD effects should not be essentially different from those 



in the SM case |f| ||. 
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of the inappropriateness of this description ensues from the fact that one expects a large 
K-factor increasing the cross-section, whereas the use of the running masses would imply 
a negative correction on the cross-section. For this reason we prefer to use the pole masses 
in our calculation, and fully parameterize our ignorance of the QCD corrections by means 
of a i^^ CD -factor whose precise value will be easily incorporated once it will be known in 
the future. 

After some digression let us come back to the general description of the cross-section 
computation for the process (|2,1|). Once a PDF for 6-quarks is used, there is some amount 



of overlap between bg- and gg- initiated amplitudes, which has to be removed 4£j. The 
overlap arises because the 6-density in the bg amplitude receives contributions from gluon 
splitting which was already counted in the gg amplitude (see e.g. the last diagram in 
Fig. ||a), so we have to avoid double counting by the subtracting of the gluon splitting 
term. The net partonic cross-section from the bg- and ^-initiated subprocesses is 

v(h + 99^ H+i+ X) net =a(gb H+t) + a{gg H + tb) 

-a{g -^bb®gb^ H + i). (2.7) 



Recalling that the 6-density evolution equation is given by the standard formula [57| 

rf W d*,Q 2 ) = ^p- |^{P^)/ 6/ft (pQ 2 )+^W/ fl/ft (^ 2 )} (2. 



d 

d\nQ 2 



where fi/h( x j Q 2 ) is the PDF for the parton i carrying a momentum fraction x at scale Q 2 
to be found in the hadron h. The subtraction term can be approximated by integrating 
the gluon-splitting part g — > b of this equation at the leading order in as{Q 2 )- Therefore 
the modified 6-density is found to be 



a s (Q 2 ) ln (Q 2 \ f 1 dz \z 2 + (1 - z^ 2 
2 7r \m 



f g/h [-,Q 2 ) ■ (2.9) 



where Pb^ g (z) = (1/2) [z 2 + (1 — z) 2 ] is the gluon splitting function. Thus the modified 
6-density contains the collinear logarithm and splitting function Pft<_ g convoluted with the 
g-density. Accordingly the expression for the subtraction term in ( |2.7D is 

Subt.term(AB -» H+tb) = 

^ I kAto.<nM^<tw,r \ (2 . 10) 

\ +f s /A(?i, Q )/s/b(»2, Q )v(gb — tH+) J 1 ' 

The subtraction term g — > bb®bg — > tH + is not negligible at all as it involves a leading 
log. We shall see explicitly its numerical significance in the next section. 

The hadronic cross-section ( |2.1| ) at LO in QCD is obtained by convoluting in the usual 
way the tree-level partonic cross-sections with the partonic densities of the two colliding 
hadrons: specifically a(qq — * H + tb) must be functionally convoluted with the light quark 
densities, also the first two terms on the right-hand-side (RHS) of (|2.7| ) are convoluted 
with the g and b -densities, and finally one sums these results and subtracts the term 
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( 2.10 ). It is seen that all these contributions to the cross-section are effectively of leading 
order a s . Moreover, from the previous considerations we will apply in the end an overall 
QCD i'T-factor to account for the NLO (or higher) effects. For definiteness we will take 
^QCD ^ 5 whenever a QCD K-factor is invoked. Once the exact value of K will be 
known, it will be easy to rescale our plots to take into account the effect of the (process- 
dependent) gluon loops. 

Also, for the study of the various differential distributions one has to properly combine 
the 2 — ► 2 and 2 — ► 3 processes in order to reproduce not only the total cross-section but 
also the correct event kinematics. The point is that we know the total amount of double 
counting but not a priori which part of this value should be subtracted from the H + t 



process and which part from the H + ib one. We apply here the method proposed in [58] 
for the analogous process of the single top quark production. According to this method we 
use the cut on the transverse momenta of the 6-quark associated with charged Higgs boson 
production to separate and recombine gg- and frg-initiated processes. We will discuss this 
in detail in the next section. As a result we obtain the correct combination of kinematical 
distributions. To the best of our knowledge this has not yet been done in charged Higgs 
boson study searches in the literature. 

Let us next briefly comment on the relevant MSSM quantum effects for the pro- 
cess ( |2.l[) . A more detailed discussion is given in Section ||. Amplitudes that are propor- 
tional to the bottom quark Yukawa coupling receive supersymmetric quantum corrections 
that can be very sizeable for large tan/3 values |33|]. This feature has been exploited for phe- 
nomenological applications to the physics of the on-shell tbH + vertex in Ref. However, 
it can also be important for the corresponding off-shell vertex involved in the production 



process (j2.1| ), as first pointed out in |45j where a first estimation was attempted. For 
tan/3 > ym^/m^ process fl2.1[ ) is dominated by the piece proportional to the bottom quark 
Yukawa coupling. Moreover, corrections of order (a/Air) n tan n /3, where a = as, aw, ap- 
pear at n-th order of perturbation theory [K^] . They can all be summed up in the definition 
of the renormalized bottom Yukawa coupling, 

h= — t— r ► — t— r tan/3 (2.11) 

v\ 1 + Am b v 1 + Am b 

where v\ = v cos (3 and so the last expression is valid only for tan/3 S> 1, showing the 
prominence of the large tan/3 region for our process. The quantity Am b in eq.( 2.11| ) is 



driven by both strong (SUSY-QCD) and electroweak (SUSY-EW) supersymmetric effects 
(depending on whether a = as or a = aw) and it increases linearly with tan/3 (Cf. 
Section ||), therefore reaching values that can be of order 1. As a result any realistic 
analysis of the reach of process fl2.1| ) in H searches clearly demands for the appropriate 
inclusion of these dominant SUSY radiative corrections. 6 Incidentally, let us recall that at 
present large tan/3 scenarios, such as those derived from supersymmetric SO(10) models 
with unification of the top and bottom Yukawa couplings at high energies [j33| , |59|| , have 
become more and more appealing since LEP searches for a light neutral Higgs boson, h, 



6 For the impact of this type of effects on the determination of the exclusion plot in (M H + , tan 0) space 
at the Tevatron, see Il7|. 
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started to exclude the low-tan/3 region of the MSSM parameter space. The latest analyses 
rule out the MSSM for tan/3 in the range 0.5 < tan/3 < 2.4, even with maximal stop 
mixing |6C]. 

We have already argued above that the charged Higgs boson production process ( |2.1| ) 
is expected to have large, and positive, QCD corrections. This circumstance, together with 
the additional cooperation from the SUSY corrections themselves, which can also be large 
and positive, could be crucial to generate a substantial "effective i^-factor" in the MSSM 
context : 

^MSSM _ ^QCD ^SUSY _ ^ .12) 

For K MSSM > 2 (e.g. through K® CD = 1.4 and if SUSY = 1.5 ) the opportunity to chase a 
MSSM Higgs boson already at the Tevatron would be open, while a more detailed analysis 
could be performed at the LHC (see Section Q). 

Concerning the method employed to compute the squared matrix elements, we have 
made intensive use of the CompHEP package |6l| , for both the signal and background 
processes. Although CompHEP is only able (in principle) to deal with tree-level calcula- 
tions, we have managed to add the supersymmetric corrections to the H + ib vertex and 
fermion propagators and we have assessed the relevance of the off-shell contributions. The 
importance and significance of the various corrections is explained in Section || where its 
numerical impact is evaluated for the signal fl2.1| ). In the meanwhile, since the background 
processes are completely insensitive to the leading SUSY effects of the type we have men- 
tioned, in the next section we present a detailed signal versus background analysis in which 
the cross-sections for all processes are computed without supersymmetric corrections. 

3. Signal and background study 
3.1 Preview 

We are interested in the search for a heavy charged Higgs boson with a mass larger than 
the top-quark mass. The case where the t — > bH + decay is allowed has already been well 



investigated for the upgraded Tevatron [g, |4j, |62j, |63J . We focus on the ttbb signal signature 



corresponding to the H + — > tb Higgs boson decay channel with the largest branching ratio. 
First we check the potential of the Tevatron collider for the heavy Higgs boson search, then 
we consider the LHC, whose detection region for M H + > mt will obviously be enlarged Q. 



Channels (|2.3|) and (^J) have been studied in 1 35 , 39] for the triple-6-tagging. Recently, 
channel ( |2,3j ) has been considered for the four-6-tagging case The four-6-tagging search 
improves the signal/background ratio but at the cost of signal rate and significance. The 
study of channel ( |2.4j ) has also been extended recently (for which both cross-sections (|2.4| ) 
and (2.3) have been combined) for the triple-6-tagging case [34|. One should note the 



additional channels that have been suggested recently for the charged Higgs boson search. 
In particular, the H^W^ [25] and H^H^ |2^] production modes followed by the decay 



— > hW^ (where h = h°,H°,A is any of the neutral MSSM Higgs boson allowed by 
phase space). However, none of these modes is favoured at high tan/3 and therefore will 
not be considered in our study. 
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In this paper we consider the dominant ttbb signature for the combined signal pro- 
cess (^j]) m detail and at a more realistic level than it was done previously. We focus 
here on the triple-6-tagging case, which gives the best possibility to measure the signal 
cross-section. This is the crucial point, especially for the Tevatron, where the production 
rate is too small to give any viable signal in the case of the four-6-tagging. In the case 
of the LHC, a triple-6-tagging study allows the signal cross-section to be measured more 
precisely (as we shall show below), even though the signal/background ratio can be better 
for the four-6-tagging case. Let us stress the following points before presenting the results 
of our analysis: 

• To study the tibb signature, it is important to take into account the effects of string 
fragmentation, initial- and final-state radiation. 



• As warned before, the processes (2.3) and (|2.4j ) should be properly combined. This 
is important because in the high p h T region those processes are qualitatively different 
(contrary to the low-p^ region). Of those events produced via gg — ► H + tb, only a 
25(40)% survive ap\ > 30GeV cut at the Tevatron (LHC). However one should note 
that the efficiency of 6-tagging is py-dependent. After the 6-tagging, the percentage 
of gg- initiated events under the p\ > 30GeV cut roughly doubles. That is why 
PT-dependent 6-tagging efficiency was used in our study. 

• We present signal and background cross-sections and give the final results in terms 
of the signal and background efficiencies. Based on this information one can derive 
the reach of the Tevatron and LHC for a given integrated luminosity. 

• Another aspect of this work is the determination of the accuracy of the signal cross- 
section measurement, which is the crucial point for the measurement of the H + ib 
coupling. This is very important in order to understand how the signal cross-section 
is affected by SUSY corrections in various regions of the parameter space. Certainly 
this is the most ambitious goal of the present study. 

3.2 Signal and background rates 

For the signal and background calculation, we use rut = 175 GeV, m& = 4.6 GeV, tan/3 = 



50, and the CTEQ4L set of PDFs [65]. Here mt, mj, refer to the quark pole masses. Unless 
stated otherwise, we assume K® CD = 1 to compare signal and background rates. Tables || 
and H present the tree-level total signal rates as well as the rates for the subprocesses and 
the subtraction term for the 2 TeV Tevatron collider and 14 TeV LHC collider respectively. 
The total signal rates are illustrated also in Fig. [2[ Note that the computation using only 
the bg — > iH + channel as done in |3?]] overestimates the cross-section by ~ 25% at the 
LHC. The subtraction term g —* bb ®bg — > iH + is indeed not negligible as it involves a 
leading log, see eqs.( 2~l)| , 2.10| ). 



The total cross-section is given only for the H + production; the inclusion of the H 
channel is just given by twice the result displayed in tables and figures. In Fig. | one 
can see that the signal rate is roughly 3 orders of magnitude higher at the LHC than at 
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M A (M H +) 


bg -> tff+ (fb) 


gg -> #+i6 (fb) 


Subt.term (fb) 


qq — ► il+ib (fb) 


Total (fb) 


200 (215) 


3.20 


1.02 


1.86 


3.86 


6.22 


250 (263) 


1.37 


0.409 


0.758 


1.05 


2.07 


300 (310) 


0.587 


0.166 


0.311 


0.353 


0.795 


350 (359) 


0.253 


0.0688 


0.128 


0.125 


0.319 


400 (408) 


0.110 


0.0280 


0.0531 


0.0466 


0.131 


500 (506) 


0.0209 


0.00470 


0.00914 


0.00689 


0.0234 


600 (605) 


0.00398 


0.000778 


0.00154 


0.00105 


0.00428 



Table 1: Tree- level signal rates at the Tevatron (2TeV) for various subprocesses and at fixed 
tan/3 = 50, if QCD = 1. 



M A (M H+ ) 


bg -> tH+ (pb) 


gg -> H+tb (pb) 


Subt.term (pb) 


qq -► H+ib (pb) 


Total(pb) 


200 (215) 


5.55 


3.03 


4.22 


0.101 


4.46 


300 (310) 


2.53 


1.30 


1.83 


0.0223 


2.02 


400 (408) 


1.22 


0.594 


0.847 


0.00742 


0.973 


500 (506) 


0.625 


0.294 


0.422 


0.00296 


0.500 


600 (605) 


0.344 


0.158 


0.222 


0.00133 


0.281 


700 (704) 


0.194 


0.0873 


0.123 


0.000641 


0.159 


800 (804) 


0.114 


0.0498 


0.0705 


0.000328 


0.0938 



Table 2: As in Table | but for the LHC (14TeV). 



the Tevatron. Also, the production rate at the LHC drops with the mass slower than 
it does at the Tevatron. For tan/3 = 50, Mjj+ = 215(359) GeV and 25 fb -1 integrated 
luminosity — 174(8) H + bosons are produced at the Tevatron. For the same value of 
tan/?, M H + = 200(600) GeV and an integrated luminosity of 100 fb -1 — 530k(32k) H + 
bosons are expected at the LHC. The dependence of the signal cross-section with tan f3 
can be appreciated in Fig. [| for both the Tevatron (at fixed M H + = 250 GeV) and at the 
LHC {Mjj+ = 500 GeV). It is obvious that the signal remains negligible in the (allowed) 
low tan P segment, say 2 < tan P < 20, as compared to its value for tan P > 30. 

We have checked the uncertainty of the signal due to the choice of PDF sets. We 
have compared the results in Tables |l] and ^ with the ones obtained with the MRST 
(central gluon) PDFs [66]. The results show a large deviation for some of the individual 
sub-channels (up to ~ 50% deviation), but they are compensated in the sum, leaving a 
5 — 10% uncertainty on the total cross-section. For example, for Mjj+ = 215 GeV (Ma = 
200 GeV) we obtain the following set of cross-sections: (bg, gg, Subt.term, qq, Total) = 
(3.80,1.63,2.95,4.52,7.01) fb at the Tevatron and (6.54,4.2,5.68,0.136,5.19) pb at the 
LHC. These values correspond to an uncertainty in the cross-section of (17, 46, 45, 15, 11)% 
at the Tevatron and (8,16,15,15,7.5)% at the LHC. The substraction procedure plays a 
key role in reducing the uncertainty, since there is a compensation between the variation 
of the gg and bg channels and the substraction term. This reduction is fully effective at 
the LHC. At the Tevatron the compensation between the gluonic channels also takes place, 
but the uncertainty is increased by that of the gg-channel, which plays a key role at that 
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Figure 2: Signal rates as a function of the charged Higgs boson mass at a) the Tevatron and b) 
the LHC (since the gluon initiated processes (2.3), (2.4) dominate at the LHC, the solid and dashed 
lines are indistinguishable in this plot), for tan/3 = 50, JfQ CD = 1. 




Figure 3: Signal rates as a function of tan/3 for a) the Tevatron with M H + = 250 GeV, and b) 
the LHC with M H + = 500 GeV, if QCD = 1. 



machine. The cross-sections obtained using the MRST PDFs are usually larger than the 
ones of the CTEQ4L. For M H + > 300 GeV at the Tevatron they become smaller, and for 
M H + = 500 GeV the MRST prediction is a 25% smaller than that of the CTEQ4L. 

We have also checked the uncertainty in the signal due to the choice of renormalization 
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and factorization scales (assumed both equal to /xr). We have analyzed the dependence 
on inside the interval M H +/2 < (j,r < 2M H +. It is obvious that a larger value of 
fiR decreases the signal cross-section. Again, individual sub-channels show a stronger 
dependence than the total cross-section. The relative uncertainty is very weakly dependent 
on the charged Higgs boson mass. For the total cross-section it is around ~ 28% at the 
Tevatron and ~ 18% at the LHC. For M H + = 215 GeV (M A = 200 GeV) we find the 
following sets of cross-sections for [ir = M H +/2: (4.07,1.73,2.52,5.63,8.91) fb at the 
Tevatron; and (5.91, 4.02, 4.66, 0.136, 5.41) pb at the LHC. For (j, R = 2M H+ the results are: 
(2.68,0.678,1.44,2.89,4.81) fb at the Tevatron and (5.38,2.35,3.87,0.0895,3.96) pb at the 
LHC. The background channels in Table || also present uncertainties due to the choice of 
PDFs and /j,r. The central value of the ^.r scale for the signal processes has been chosen 
equal to M H +, whereas that of the background processes to (2m t ). 

As we have already mentioned we will focus on the tibb signature. We consider the 
case where one top decays hadronically and the other leptonically (including only electron 
and muon decay channels) in order to reduce the combinatorics when both top quarks are 
reconstructed. The branching ratio of tibb — > bbbbt^vqq' is 2/9 x 2/3 x 2 = 8/27. 

In order to decide whether a charged Higgs boson cross-section leads to a detectable 
signal, we have to compute the background rate. Since the miss-tagging probability of light 



quark and gluon jets is expected to be ^ 1% |38, 69], the main backgrounds leading to the 



same tibb signature and their respective cross-sections are those shown in Table y. 

The above notations qq, gq, gb for background processes assume that we have summed 
over initial-state quarks and anti-quarks. One should notice that the double counting, 
as in the signal case, appears also when one sums the gg — > tibb and gb —* tib processes. 
Therefore one should subtract the overlap between them, which we denoted by Subtr. term 
- Cf. eq.(|2l0|). 



In considering our final state signature tibb we should point out that the partial width 
of the decay mode H + — > tb can be itself sensitive to important SUSY radiative correc- 
tions [ p.2jj . However, in the present instance it is only the branching ratio of this process 
that enters the calculation. At high tan/3, the only relevant mode other than H + — > tb is 
H + — > t + v t and the latter is of order 10% at most. Therefore BR(H + — > tb) is not too 
sensitive to SUSY effects. By the same token BR(H + — > t + u t ), though smaller, can be 



quite sensitive [12] and so with sufficient statistics it could be used as an additional test of 
the underlying SUSY physics. We factorize these corrections from the production process 
itself, and we will take them into account only in the combination of the signal/background 



analysis and the radiative corrections in Section 4.3. 
3.3 Simulation details 

To perform a realistic signal and background event simulation we complied to the following 
procedure. The matrix element for the complete set of signal and background processes 
has been calculated using the CompHEP package |HJ. The next step was the parton- 
level event simulation, also with the help of CompHEP. Then we automatically linked the 
parton-level events from CompHEP to the PYTHIA6.1 Monte Carlo generator |7C|] , using 
the CompHEP-PYTHIA interface jfl|. 
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(a) pp(pp) -► ttb(b) 



Tevatron 



a{qq — > ttoo) = 


6.62 lb 






0.676 fb 




a {go — ► tto) — 


l.ZZ ID 




Subtr. term : 


0.72 fb 




Total : 


7.80 fb 










(T{QQ — ► ttOO) — 


U.ZDD pu 




<7\99 ~^ ttb0 ) — 


D.UU pu 




(J [go — > tto) — 


4.00 pu 




Sm 1 n t" Y T f^VYH 
UUUtl > UCJ-111 


2.1 pb 




Total : 


8.50 pb 




(b) pp(pp) -> ttqg 


when the light quark or gluon 


are misidentified as a b-jet 


Tevatron 


LHC 




a(qq -> gtt) = 


1890 fb a(qq -> gtt) = 


21 pb 


a{gq -> gtt) = 


193 fb a(gq -> qtt) = 


122 pb 




262 fb a(gg gtt) = 


371 pb 


Total : 


2345 fb Total : 


514 pb 



Table 3: The main background processes for the charged Higgs boson production at the Tevatron 
and LHC. For the ttbb and ttqg processes we have applied the jet separation cut Ajj > 0.5(A^ = 
y/A6 2 + A0 2 ) and the cut jP T > 10 GeV (p 3 T > 20 GeV) at the Tevatron (LHC). For the ttj process 
the cut jP T > 10 GeV (p 3 T > 20 GeV) was applied at the Tevatron (LHC). 



Therefore we took into account the effects of the final-state radiation, hadronization 
and string-jet fragmentation using PYTHIA tools. The following resolutions were used 
for the jet and electron energy smearing: AE had /E = 0.8/VE~ and AE ele /E = 0.2/y/E. 
In our analysis we used the cone algorithm for the jet reconstruction with a cone size 
AR = \J Alp 2 + A-rf = 0.7. The choice of this jet-cone value is related to the crucial role 
of the final-state radiation (FSR), which strongly smears the shape of the reconstructed 
charged Higgs boson mass. We have checked that the value of 0.7 minimizes the FSR 
effects. 

The minimum Et threshold for a cell to be considered as a jet initiator was chosen to 
be 5 GeV (10 GeV) for the Tevatron (LHC), while the minimum threshold for a collection 
of cells to be accepted as a jet was chosen as 10 GeV and 20 GeV, respectively for the 
Tevatron and the LHC. 

As already noted, we require three b-jets to be tagged. A realistic description of the 
b-tagging efficiency is therefore very important. In the case of the Tevatron, we use the 
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projected 6-tagging efficiency of the upgraded D0 detector |]6q] : 

£6 = .57^ to „h(^), (3.1) 
For the LHC, we parameterize numerical results from the CMS collaboration [69]: 



0.6, for p T > 100 GeV 

0.1 +p T /(200GeV), for 40 GeV < p T < 100 GeV (3.2) 
1.5pr/(100GeV) - 0.3, for 25 GeV < p T < 40 GeV 



We assume that 6-jets can be tagged only for pseudorapidity |%| < 2 by both Tevatron 
and LHC experiments. 

3.4 Combining the bg — > H + i and gg — > H + tb processes 



As we mentioned above, we apply the recipe of Ref. [58| to combine bg — > H + i and 
gg — > H + ib in order to get the correct overall distributions. 

The tree level 2 — > 2 and 2^3 processes reproduce correctly the px distribution of the 
associated 6-quark (p T ) only in certain parameter regions. The {p\) distribution of 2 — > 2 
(for which the 6-quark comes from initial state radiation simulated by PYTHIA) is correct 
only for small p T values, because the gluon splitting function is not able to reproduce high 
p T region correctly. Contrary, the 2^3 processes reproduces correctly the distribution at 
high p h T since it includes the complete set of respective diagrams, but it fails to reproduce 
the correct p T in the low Pt region where one should take care of the resummation of large 
values of log[M H+ /M b ]. 

We have compared various kinematical distributions of the 2^3 process pp(gg) — > 
H + ib and 2 — > 2 process pp(bg) — ► H + i+b sp nt- We have found a proper matching between 
the resummed contribution in the collinear region for the 6-quark and the complete tree- 
level contribution in the hard region. Figure |I| shows transverse momenta and rapidity 
distributions for the final state H + , t-quark and 6-quark for the processes at the Tevatron. 
For the LHC the distributions are qualitatively the same. As expected the difference is 
clear in the 6-quark distributions. For the pp(bg) — > H + i + b sp u t process, the 6-quark is 
softer and less central than that for the pp(gg) — ► H + tb process. At the same time it is 
shown that the i? + -boson and i-quark distributions are nearly the same. 

We use the method of matching collinear and hard kinematical regions based on the 
kinematical p\ separation of the pp(bg) — > H + t + b sp ut and pp(gg) — > H + tb processes in 
the regions p\ < p^ 1 and p\ > p^, respectively. We search for the value of p™ f in order 
to satisfy two requirements, namely: 

1) the common rate of pp(bg) — > H + t + b sp u t with p\ < p™* and pp(gg) — > H + tb with p h T > 
p™ 1 gives the combined total rate computed in the previous section; in other words, one 
can normalize a rate in a collinear region on the Ototai ~ &\pp{99 + QQ) ~^ H + ib, p h T > p™*]; 

2) the overall p h T distribution should be smooth. 

The result is illustrated in Fig. ||, where we show several variants of the combination of 
these two processes for various values of p™ 1 . We found that the optimal p™ 1 providing a 
smooth sewing for these two processes at the Tevatron and LHC is equal to about 25 GeV. 
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This value gives physically reasonable answers on the main questions of this section: 

a) in which kinematical regions should the pp(bg) — > H + i + b sp ut and pp(gg) — ► 
processes be considered and how one should properly simulate them; 

b) how the subtraction term is distributed between pp(bg) — > H + t + b sp ut and pp(gg) — > 
H + tb processes, and what part of double counting should be subtracted from each subpro- 
cess. One should notice that in some particular cases, like the one chosen for illustration 
in Fig. [|, there is practically no difference in choosing p^ 1 in the range of 25 — 40 GeV. 
This can be seen from Fig. || as well as confirmed by our numerical results for the final 
efficiencies. 

We conclude that the method of combining the pj, distribution of H + t + b sp m and the 
complete tree-level H + ib process allows us to find the physically motivated pr cut on the 
6-quark, which allows us to treat together those processes and simulate them in different 
kinematical regions of p b T . Namely, we generate pp(bg) — > tH + + b sp nt events using the 
PYTHIA generator, and use them in the kinematical region p h T < p™* with the weight, 
corresponding to the a to tai — &[pp{gg) — * H + tb, p h T > p™*] cross-section; in the region 
p\ > p™ 1 , instead, we use pp(gg) — ► H + ib events. 

3.5 Kinematical analysis 

As we mentioned, the bbb(b)+£ ± +2jets+p'!p iss signature of signal and background processes 
is the subject of this study. One should reconstruct the ttb(b) state from this signature 
and then, the charged Higgs boson mass from all possible combinations of i6-invariant 
masses. At the same time one should work out an efficient set of kinematical cuts for the 
background suppression. 

For reconstructing the tib(b) final state, we follow the procedure described below. 

• We reconstruct the W-boson mass from lepton and neutrino momenta: = 
(pe + Pv) 2 ■ The basic cuts for the lepton (electron or muon) has been chosen as 
follows: 

]4>15GeV, |%| < 2.5, p^ ss >15GeV. (3.3) 

To find the p z of neutrino and therefore the neutrino 4-momentum, one should solve 
the quadratic equation = M]y, which can have two solutions. We reject events 
if this equation has no solutions, while in case it has two solutions we keep both of 
them. 

• We reconstruct the mass of the second W-boson (M^) : we keep all dijet combina- 
tions for the moment: the effects of the initial- and final-state radiation and 6-miss- 
tagging, the number of light jets is almost always larger than two. The following 
basic cuts for the jets were chosen: 

p^ b > 20 (30) GeV [Tevatron(LHC)] , \r]j\ < 3, \rj b \ < 2. (3.4) 

• Then we form and keep all mt 1 = Mwib an d mt 2 = My/2b combinations for the first 
and second top-quarks. 
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Figure 4: Transverse momenta and rapidity distributions for the final-state particles of the pro- 
cesses pp(bg) — ► H + i + b S piit (dashed line), a 2 — ► 2 process with an additional 6-quark from 
initial-state radiation, and the 2^3 pp(gg) — ► process (solid line) at the Tevatron for 

M H+ = 300 GeV. 



• In the final step we form the \ function 



X = - M w y + (M^l - M W Y + (m tl - m t y + K - m t ) 2 (3.5) 
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Figure 5: Transverse momenta distribution of the 6-quark sewn for various values of p™ 4 for the 
pp(bg) — ► H + t + b spin and pp{gg) — > H + tb processes at the Tevatron for M H + = 300 GeV. 

for all combinations of 6-jets, jets, lepton and neutrino and choose the combination 
giving the smallest (best) value of the \ function. The optimized value of the cut on 
the x function was found to be 

X < Xmax = 100 GeV at the Tevatron, % < Xmax = 60 GeV at the LHC . (3.6) 

It should be noted that the reconstruction thus made is independent of the order in which 
Mp^, My^2-> m ti > fnt2 were formed. This leads to a better signal efficiency, a better probabil- 
ity of a correct reconstruction and a better control of the efficiency through the only relevant 
parameter, Xmax- At the Tevatron, the cut on x is quite relaxed by the lack of statistics, 
while at the LHC this cut could be tightened further, gaining both signal/background 
(S/B) ratio and significance (S/y/B). 
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After the reconstruction of the ttbb state one should reconstruct the charged Higgs 
boson mass for the signal and the continuous tb mass for the background. We assume that 
the 6-jet with the highest px in tbt(b) signature comes from the H + decay. The probability 
of that being correct is directly related to the value of M H + : for M H + ~ 200 GeV it is only 
about 50% while for M H + ~ 300 GeV the 6-jet coming from the H + decay has the highest 
Pt in already 75% of the cases. Since we chose just one 6-jet, there will be two ifr-invariant 
mass combinations with two top-quarks that will enter the rritb invariant mass plot. 



Like in Ref. [34|, we found that p\ from H + — > tb decay is a good variable for the 
separation of signal from background. However, instead of the fixed p h T cut we apply here 
the -dependent p\ cut: 

p b T > [M H + /5 - 15] GeV , (3.7) 

since the peak of this 6-jet distribution is just proportional to M H + . This M H + dependence 
of the p\ cut allows us to increase the efficiency of the kinematical cut and selection for the 
signal. The cut depends on the Higgs mass and should be understood as the cut chosen 
with the respect to the mean value of the mass window where we are looking for the Higgs 
signal. After the fitting of the Higgs mass using the 'rough' pj,, cut value, one can use the 
fitted Higgs mass as the input for p b T {M H +). 

The window around the selected values of M H + was also chosen M^+-dependent: 



\m tb - M H+ 1 < 5^M H+ . (3.8) 

One could think of an additional set of kinematical cuts, such as the correlation angle 
between the Higgs boson direction and its decay products (in the Higgs boson rest frame), 
the angle between top and bottom from Higgs boson decay (since they tend to be more 
back-to-back). It turns out that already at the PYTHIA simulation level the difference 
in those distributions between signal and background is quite blurred. The application of 
the respective cuts would lead to some increase of the S/B ratio, but at the same time 
to a definite decrease of the significance and of the accuracy of the signal cross-section 
measurement. 

After all cuts are set up we are ready to present signal and background efficiencies, 
S/B ratio and significance. The results are summarized in Tables || and |5| for the Tevatron 
and the LHC respectively. We present there the number of signal (^^)-( |2^ ) [S] events 
and of background (Table |||) [B] events, as well as their ratio and significance after the 



reconstruction procedure, 6-tagging (3.1)- (|3.2|) and set of cuts (p^3|)- (|3.&j ). Numbers are 
given for tan/3 = 50 for an integrated luminosity of L = J Cdt = 25 fb _1 at the Tevatron and 
an integrated luminosity of 100 fb _1 at the LHC. The number of signal events corresponds 
the tree- level cross-section (K QCD = 1). Tables §§] also present the various efficiencies (e) of 
the cuts ( |3 . 3| ) - ( J378D and reconstruction (including 6-tagging) for the signal and backgrounds. 
The last columns of these tables give the 95% C.L. (resp. 5er) discovery limits of the total 
signal cross-section in femtobarns (resp. picobarns) for a 25 fb _1 (resp. 100 fb _1 ) of total 
integrated luminosity for the given efficiencies at the Tevatron (resp. LHC). 

As an example, Figure ^ shows the reconstructed tb invariant-mass distribution for 
signal and background events at the Tevatron and the LHC. Owing to the much higher 
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Figure 6: Reconstructed tb invariant-mass distribution for signal, background and signal plus 
background events, for a) Mh+ = 263 GeV (-M40 = 250 GeV), at the Tevatron and b) Mjj+ = 
506 GeV (M A a = 500 GeV), at the LHC, if QCD = 1. 



M H+ (GeV) 
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S/B 


s/Vb 
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95% C.L. (fb) 


215 
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14.0 


0.70 
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0.060 


7.50 


263 


3.5 


7.7 


0.46 


1.27 


7.70 


3.10 


0.034 


5.39 


310 


1.6 


7.7 


0.21 


0.58 


8.90 


3.10 


0.034 


4.66 


359 


0.7 


6.4 


0.10 


0.26 


9.20 


2.70 


0.028 


4.29 



Table 4: Number of signal ( 2.2 )-( 2.4 ) [S] and background -Table fs- events [B], their ratio and 



significance after reconstruction procedure, 6-tagging (3.1)-(3.2) and set of cuts (3.3)-(3.8) at the 
Tevatron. The last column gives the 95% C.L. discovery limit (~ 2 a) on the total signal cross- 
section in femtobarns for 25 fb -1 of total integrated luminosity and given efficiencies. See text for 
details. K^ CD = 1. 



signal rate, the significance of the signal at the LHC is considerably higher than at the 
Tevatron. In addition, the difference in shape of the signal and background distributions 
could be clear only for M H + > 400 GeV, which will be accessible only at the LHC. 

In spite of the fact that the S/B ratio is quite high at the Tevatron for M H + = 215 GeV, 
the signal rate itself is small, even for tan/3 = 50. That is why we try to keep the signal 
efficiency as high as possible by relaxing the acceptance cut on p!^ and the cut on \- 
Table |] shows that at the Tevatron there is no 5 a discovery limit available for Higgs boson 
masses above 200 GeV, although there is a narrow window (175 — 195) GeV for Higgs boson 
masses just above the kinematical limit for the top quark decay into a charged Higgs boson. 
In general, however, at the Tevatron one can only hope to exclude charged Higgs boson 
masses up to the certain value defined by the efficiencies, the value of tan (3 (and the other 
MSSM parameters beyond the tree level) and the total integrated luminosity. Even the 
exclusion region is apparently small. The signal significance is only 3.4 for the charged 
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iV%+(GeV) 


S 
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215 
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0.57 
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0.005 


1.02 
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31U 
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rnnn 
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0.41 


31.1 


Z.UO 


L.ZO 


U.U18 


U.32 
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1560 
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0.25 


19.8 


2.70 


1.50 


0.016 


0.25 
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773 


3770 


0.21 


12.6 


2.60 


0.97 


0.009 


0.20 


605 


433 


2070 


0.21 


9.5 


2.60 


0.56 


0.004 


0.15 


704 


217 


1170 


0.19 


6.3 


2.30 


0.33 


0.002 


0.13 


804 


117 


666 


0.18 


4.5 


2.10 


0.19 


0.001 


0.10 



Table 5: Number of signal ( 2.2 )-(pT^) [S] and background -Table ||- events [B], their ratio and 
significance after reconstruction procedure, 6-tagging ( 3^1 )-(3 1 2 ) and set of cuts (^3)- (3J3) at the 
LHC. The last column gives the 5a discovery limit on the total signal cross-section in picobarns for 
100 fb _1 of total integrated luminosity and given efficiencies. See text for details. K^ CD = 1. 




Figure 7: The total cross-section a(pp(pp) — > H + ib + X) for fixed tan/3 = 50 including a QCD 
factor ifQ CD = 1 (solid) and K^ CB = 1.5 (dotted): a) at the Tevatron, dashed and dot-dashed 
lines correspond to the cross-section values necessary for the 2a (95% C.L.) exclusion and the 3a 
evidence respectively, at the integrated luminosity L = 25 fb~ ; b) at the LHC, dot-dashed and 
dashed lines corresponding to the cross-section necessary for the 5a discovery at the integrated 
luminosity L = 30 fb _1 and 100 fb _1 respectively; the long-dashed line at the bottom corresponds 
to the 95% C.L. limit at L = 100 fb -1 . 

Hiffes boson with mass 215 GeV, tan/3 = 50 and K-factor K® CB = 1.5. For this value of 
tan/? the 95% C.L. limit on the charged Higgs boson mass is 240 GeV. We use the Poisson 
statistics when putting this limit at the Tevatron, since the number of events is small for 
both the signal and the background. On the other hand, for i^Q CD < 1 it is not possible 
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to obtain any limits using this channel, as it is obvious from the comparison of Tables | 
and |. 

At the LHC, the situation is much better, with higher luminosity (30 — 300 fb _1 ) and 
a signal rate higher by 3 orders of magnitude with respect to the Tevatron. Here we can 
allow ourselves to use the tighter kinematical cuts ( |3.3| )-(3,6). We try to reach the best 



significance ratio S/vB while keeping the S/B ratio above the 10% level. Based on the 
efficiencies of Table |5j we find that a charged Higgs boson can be discovered at the 5a level, 
with mass up to about 900 GeV for tan/3 = 50 and i^Q CD = 1.5. 

It is important also to understand the role of the i^-factor. In Fig. [?] we present the 
cross-section for the pp(pp) — > H + tb + X process for K QCB = 1.5 and if QCD = 1 (dotted 
and solid lined, respectively), and the cross-section for necessary for the discovery/exclusion 
of a charged Higgs boson. In Fig. |7|a we show the 95% C.L. limit (dashed) and the 3a 
evidence cross-section at the Tevatron for an integrated luminosity of 25 fb~ . Fig. [?]b shows 
the cross-section for the 5a discovery (i.e. 5 yB) at the LHC for an integrated luminosity 
of L = 30 fb^ 1 (corresponding to 3 years of low-luminosity regime) and 100 fb _1 (1 year 
high-luminosity regime). Notice that the discovery curves also evolve with M H + because 
the efficiencies are M H + -dependent (see Tables || and |5|). Even at low luminosity, the 5 a 
discovery limit for the LHC is 600 GeV for K® CD = 1, 700 GeV and for K® CB = 1.5. Of 
course, if the iCQ CD factor turns out to be below 1, the mass reach is lowered accordingly, 
e.g. 450 GeV for i^Q CD = 0.5. One can see that the value of the iT-factor is really 
important because a variation of j^Q CD from 1 to 1.5 causes a 100 GeV shift of the limit 
on the charged Higgs boson mass. For L = 100 fb" 1 the discovery limit is 790 GeV if 
JfQCD = i ; anc i 900 GeV for K QCD = 1.5 (600 GeV for K QGD = 0.5). In this figure we 
also show the 2 \fB curve, from which one obtains immediately that the charged Higgs 
boson can be excluded at the 95% C.L. up to a mass of about ITeV even for = 1 

at 100 fb -1 of integrated luminosity. In Section ||| we shall come back to this figure and 
we will show that there are regions of the MSSM parameter space where the total effective 
^mssm f ac t or m eq.( 2.12 ) could be ~ (40 — 50)% larger than the typical X-factor expected 



from pure QCD alone. 

One can assume that systematic uncertainties due to the parton distribution uncertain- 
ties as well as to the uncertainty from the higher-order QCD corrections for the background 
(or from a choice of the QCD scale for the background) could be significantly reduced by 
normalizing the MC tb distribution to the data. Under the assumption that normalization 
to the data will reduce systematic uncertainties up to the level of 7%, one can also estimate 
the accuracy of the measurement of the signal cross-section using efficiencies from Table || 
For example, for M H + = 408 GeV the accuracy of the cross-section measurement will be: 



8S _ V i^B) 2 + SB syst _ ^6210 + (0.07 X 6210) 2 ^ oq0/ 



S S 1560 

This means that H + ib coupling could be measured with an accuracy of (1/2)5S/S ~ 14%. 
In the next section we show that the SUSY corrections could be substantially higher than 
14%, therefore being potentially measurable! 
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4. The role of the SUSY corrections 



4.1 The leading effects. Theoretical discussion 

In Section |3| we have concentrated on the results for the cross-section and background for 



the charged Higgs boson production process (2.1) without including the various sources of 



potentially important SUSY corrections. However, in Section y we have already warned 
about the possibility of non-negligible SUSY effects in the signal and we have correctly 
identified the region of the parameter space where they can be optimized, namely at high 
tan (5. It is now time to study this issue in detail because, as we shall see below, the super- 
symmetric loops can be very sizeable and may play a momentous role in the production of a 
charged Higgs boson. In fact, they can not only enhance the signal versus the background, 
but they can also be used as a means to characterize the super symmetric nature of the 



charged Higgs boson potentially found in hadron colliders through the mechanism (2.1). 

Among the plethora of possible SUSY corrections we disregard virtual super symmetric 
effects on the gqq and ggg vertices and on the gluon propagators. We expect those to 
be of order (a s /47r) • (s/Mg USY ) and thus suppressed by a non-enhanced (i.e. tan/3- 
independent) MSSM form factor coming from the loop integrals. Therefore, we can neglect 
these contributions as we are only considering effects of the form (a s /47r) n • tan n at 
large tan/3 . As previously emphasized, the cross-section for the signal increases steeply 
with tan/? (see Fig. |3|) and becomes highly significant for tan/3> 30, while it is much 
smaller for tan/3 in the low interval 2 — 20 where the remaining SUSY corrections are of the 
same order or even dominant, so our approximation is well justified. Similarly, we neglect 
all those electroweak corrections in vertices and self-energies which are proportional to 
pure SU(2)l x U(1)y gauge couplings; in particular, vertices involving electroweak gauge 
bosons and those involving electroweak gauginos 7 . Furthermore, we have checked that 
vertices involving Higgs bosons exchange yield a very tiny overall contribution, due to 
automatic cancellations arranged by the underlying supersymmetry. Finally, there are the 
strong gluino-squark diagrams and the tan /3-enhanced higgsino-squark vertices implicit in 
chargino-neutralino loops. We have extracted (see below) the parts of these interactions 
which are (by far) the more relevant ones at high tan /? and confirmed that the remaining 
contributions are negligible. In practice this means that we will concentrate our analysis 
on the interval tan (5 > 20 where we can be sure that our approximation does include the 
bulk of the MSSM corrections while at the same time the cross-section of process ( [2.1| ) 
starts to be sufficiently large to consider it as an efficient mechanism for charged Higgs 
boson production (Cf. Fig. ||). To be more precise, we shall hereafter confine our study 
of process ( |2.1| ) to within the relevant interval 20 < tan/3 < 60, where the upper limit on 
tan/3 is approximately fixed by the condition of perturbativity of the Yukawa couplings. 

As a consequence we can just concentrate on the leading quantum effects. The latter 
can be conveniently described through an effective Lagrangian approach that contains 
effective couplings absorbing both the leading SUSY contributions and the known part of 



7 We have performed the full calculation in the squark and chargino-neutralino mass-eigenstate basis, 
but we have cross-checked the identification of the leading parts using also the weak-eigenstate basis, i.e. 
in terms of diagrams involving squarks, gauginos and higgsinos. 
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the QCD corrections [35]. At high tan/? the most relevant piece is the effective tbH + - 
coupling as it carries the leading part of the quantum effects. Indeed, on the one hand 
the quantity Am b in ( 2.11| ) contains the bulk of the super symmetric contributions. On 



the other by trading the bottom quark mass by the corresponding running quantity one 
can also include a sizeable part of the QCD effects, namely the (universal) renormalization 
group effects. As a result the relevant piece of the effective Lagrangian at high tan /3 
involving the tbH + vertex can be cast as follows: 

r gV tb m b {^ R ) tan /3 . 

£.= —= H^t L bR + h.c. (4.1) 

y/2M w l + Am b y ' 

where fiR is the renormalization scale for the bottom quark running mass fn b (^LR) at the 
NLO in the MS scheme. One usually assigns a value to given by the characteristic 
energy scale of the process under consideration. In our case (Cf. Fig. 1) it would be 
natural to set //r = nit + M H + . However, following the discussion in Section ^, we will use 
the (on-shell) pole quark masses in the tree-level cross-section, and will parameterize the 
QCD corrections by means of a i^^ CD -factor. 

On the SUSY side the effects encoded in Am b are related to the bottom mass coun- 
terterm Jmj in the on-shell scheme. If m b is the pole mass, then its relation with the 
corrections resummed into the bottom quark Yukawa coupling leads to the consistency 
formula 

m b + 5m b =— — ; (4.2) 

1 + Am b 

in particular Am b = — 5m b /m b at one loop. The SUSY effects on Am b can be of two types, 
SUSY-QCD and SUSY-EW. The strong part of Am b originates from finite corrections 



induced by mixed LH and RH weak-eigenstate sbottoms and gluino loops [1C, |33| (Fig. ||a): 



, x n as (Msxjsy) b 
( Ara WsusY-QCD = - °F ^ m-g M LR I{m~ bi , m-^ , m gJ 



2a 5 (M SUSY ) , , , , 

~~ *■ H ^- m g ^ tan P I(m bi ,m b2 ,rn g ) , (4.3) 

where Cp = (N^ — V)/2N C = 4/3 is a color factor (with N c = 3). Moreover we have 
as(MgusY) — 0.09, M\ R = A b — fx tan/3 and the last expression in (|4.3|) holds only for 
sufficiently large fi tan/3 3> A b . Here fi is the higgsino mass parameter in the superpo- 
tential, and should not be confused with the renormalization scale fiR. Finally, A b is the 
bottom quark trilinear coupling in the soft-SUSY-breaking part of the superpotential. As 
for the electroweak effects on Am b they stem from similar loops involving mixed LH and 
RH-stops and mixed charged higgsinos #i - H 2 and read (Fig. |b) 

( A ™f>)sUSY-Yukawa = +^ ^ ™t M^I^ , m t ~ , (J,) 

Y 

-> +~r M H tan/3 I (mi , m l \i) . (4.4) 

47T 

Here At is the top quark counterpart of A b ; we have defined Yt = ^f/47r ~ 1/47T and 
M l LR = At — fj, cot (3, and again the last expression holds for large enough tan/3 only. 
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Figure 8: Leading contributions to Am b - eq. (4.2). a) SUSY-QCD contributions (4.3); b) 



SUSY- Yukawa contributions (4.4); c) SUSY-EW gauge contributions (15). Each cross represents a 
mass-insertion; each cross with a circle represents a field mixing term in the electroweak eigenstate 
basis for squarks or gaugino-higgsinos. 



Gaugino-higgsino (W — H 2 ) mixing also contributes with tan/3 enhanced terms to Am b 
(Fig. |c): 



(Am b ] 



SUSY-gauge 



16vr 2 



fi tan f3M 2 



cos 2 9- t l(m\ , M| , /x 2 ) + sin 2 Q~ t I(m\ , M 2 2 , ^ 



1 



+ - cos z Or I (mi , Mi,^) + sin z Or I {ml , A#, 



2 -o-v -fe!'--^"- / 1 ~ -o-\-"6 2 

In the above formulae we have introduced the (positive-definite) three-point function 
I(mi,m 2 ,m 3 ) = 16 7r 2 i C (0, 0, mi, m 2 , m 3 ) 

2 2 i m ? 1 2 2 1 Tl? 1 2 2 1 

mj In — 3- + m 2 In — g + m| TO3 in 



(4.5) 



(mf — m 2 ) (m 2 — 777.3) (m\ — 777.3) 



(4.6) 



This function is well defined even when all masses are equal, which is precisely the situation 
where it attains its minimum. Thus, if we roughly assume that all masses are equal to 
Msusy, it boils down to 



I(m; = Msusy) 



1 



9 M 2 

z 1VI SUSY 



(4.7) 



Furthermore, by setting \fi\ = m~ g = Msusy and tan/3 = 50 one immediately sees that 
eq.( f4.3|) gives a large correction of ±0(50%) depending on the sign of //. Notice that in the 
on-shell renormalization scheme defined in flCfl, the SUSY-EW gauge contributions from 



eq. (4.5) are compensated for by those equivalent to Fig. ||c but with a external r-lepton 
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and sneutrinos/r-sleptons circulating in the loop. 8 However the resummation in ( |4.1| ) 



can only be done in a mass- independent renormalization scheme [35], and the diagrams 
in Fig. |8]c must be taken into account. The difference between the two renormalization 
schemes amounts to a finite shift in the definition of tan/3 which, however, is not significant 
in the interesting region tan/3 3> 1. 

Of course, the top quark Yukawa coupling receives Amt corrections which are the 
counterpart to those for the bottom quark, eqs. (|4.3| )-( p3[ ), but they are suppressed (rather 
than enhanced) by tan/3 and so they will be ignored. 

It should be emphasized that both types of SUSY contributions to Ami, can be sizeable; 
naively one could think that the SUSY-EW contributions induced by higgsino-gaugino 
exchange are suppressed with respect to the (gluino induced) SUSY-QCD contribution by 
a factor of aw /otg ~ 10 -2 , but this is not necessarily so. In fact, in regions of the parameter 



space where the 3-point functions on the RHS of (4.3) and ( |4.4D are approximately equal 9 
the two types of corrections are roughly of order as fJ- tan (5/Msusy and Y t \x tan (3/Msusy 
respectively, and since as — Yt they should be comparable. However, from a more accurate 
estimation including the color factor Cp the ratio between the SUSY-QCD and SUSY-EW 
effects yields 

' as \ I m 



Y t ) V A t ) \A t ) 1 

This ratio can be smaller than 1 if At > 3 rrig. However in the numerical analysis below we 
typically choose m g > A t (Cf. Table |6|), and this makes the SUSY-QCD effects generally 
dominant. There are indeed reasons to envision this kind of scenario. First, the gluino mass 
is expected to be rather large, of the order of 1 TeV. Second, we cannot choose At too 
large, say above ITeV, because this could disrupt the SU(3) C and U(l) em symmetries of 
the electroweak vacuum state. An approximate necessary condition for this not to happen 
is f7| 



A 2 t < 3 (Mq + My + M 2 H2 + M 2 ) . (4.9) 

Here Mqjj are the LH and RH stop SUSY-breaking mass parameters, M^ 2 (with Mfj 2 < 0, 
due to the radiative breaking of the gauge symmetry) is the soft-SUSY-breaking mass 
associated to the H2 doublet (the one that couples the top quark) and [i is again the 
supersymmetric higgsino mass parameter. In a typical scenario where rrig is larger than 
any squark mass and any electroweak mass parameter, one may assume At < Mn t jj < rrig 
(Cf. Table g). In this case the RHS of eq.(g]§) should be greater than 1 and the SUSY-QCD 
effects are dominant. Third, these effects have the remarkable property that they decouple 



very slowly with the gluino mass [1C]. Indeed, for sbottom masses of a few hundred GeV 
the gluino mass can be as high as 1 TeV and the correction is still of order (20 — 30)%. 
This is due to the helicity flip in the gluino line which generates the gluino mass term in 



See Fig. 21b and eq. (86) in the first reference of |1C 



9 Notice that the function /(mi, 7712,7713) in (4.6) is slowly varying unless there are large hierarchies 
among the masses. 
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the numerator of eq, (|4.3| ). For m g ~ 1 TeV it would be contrived to assume At > 3wij, 
and hence the SUSY-EW effects must surely have gone away in the heavy gluino regime 
while the strong supersymmetric effects can still remain. Of course the correction ( |4.5| ) 
eventually dies out with rrig, but very slowly. Notice that in the opposite limit, m g —* 0, 
the SUSY-QCD effects (4.3) also tend to zero and one is left with the subleading SUSY- 



QCD contributions from the renormalized three-point functions. We shall not focus on 
the light gluino case here, although in practice it is equivalent to assume that the bulk of 
the SUSY corrections is just given by the electroweak part (fQjj), which can still be rather 
substantial. These remarks show that the leading type of SUSY effects on our process can 
have different origin and remain significant in large portions of the MSSM parameter space. 

Ultimately the origin of the exceptional SUSY quantum effects fl4.3|) and ( [4.4D stems 
from the breaking of SUSY through dimensional soft terms like gaugino masses and trilinear 
couplings, and therefore they cannot be described - unlike the conventional QCD effects 
considered before - through the universal renormalization group running of the parameters. 
Technically, they are finite threshold effects which must be included explicitly at energies 
above Msusy, the characteristic SUSY breaking scale at low energies pq ]. 

Remarkably, the finite threshold effects do not vanish when the soft-SUSY-breaking 
parameters, together with the higgsino mass parameter fi, are scaled up simultaneously to 



arbitrarily large values. This is immediately seen from (|4.3| ) and (4^) by letting simulta- 
neously fi —* Xfi and fh —* Am in these equations, where A is an arbitrary dimensionless 
scale factor and fh is any soft-SUSY-breaking parameter. However, it is important to no- 
tice that the apparent non-decoupling behaviour exhibited by this kind of SUSY effects is 
rendered innocuous after resumming all the corrections of the form (as ^ tan f3 /Msusy) 12 
and (Y t /j, tan (3/ Msusy)™ to all orders n = 1,2,... |35|. The resummed result is indeed 
given by the effective bottom quark Yukawa coupling ( |2.11| ), which shows that there are 
no higher order disturbing effects of the form Am b n ~ (a s /4Tr) n • tan n /3 when Am b > 1 
and the expansion 

i 

1 - Am b + Am b 2 - ... (4.10) 



(1 + Am b ) 



is not possible. On the other hand for |Am. | <C 1 the expansion is of course allowed and 
the higher order powers Am b n give small corrections to the linear approximation 1 — Am b 



used in [10]. However, for (Am-J < 1 and Am b < the resummed SUSY contribution 
is seen to substantially reinforce the one loop result. For example, if Am b = —50%, the 
one-loop correction is +50% whereas the resummed result gives a +100% effect! Therefore, 
the SUSY radiative corrections can be quite large, and may remain so even for very high 
values of Msusy, but their resummation is eventually kept under control when |Am 6 | S> 1. 

The effective Lagrangian approach just outlined is also very useful for the practical cal- 
culation of the cross-sections. In fact, by appropriately modifying the MSSM Feynman rules 
using the effective vertex interactions described above we have managed to automatically 
generate the leading quantum effects from the CompHEP calculation of the cross-sections. 

Although Am b in the above equations is the only correction of order (a/47r)tan/3 
(q = as, olw) that dominates for large tan/3, we have also tested the effect from the off-shell 
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SUSY-QCD and SUSY-EW corrections to the H + tb vertex and to the fermion propagators. 
As already mentioned above, some of these off-shell vertices (in particular those involving 
chargino-neutralino exchange) carry some tan /3-enhanced Yukawa couplings that could be 
important. Nevertheless it turns out that these contributions are not comparable to the 
leading Am b ones. The point is that the effects entering Am b at high tan/3 are of the 
type as tan (3 or Y t tan (3 whereas those from the aforementioned vertices are of the type 
(g nib tan (3 /My/) 2 /47r ~ Yt, and so they are down by an extra power of tan (3. At the end 
of the day they are subleading effects on equal footing to the previously dismissed SUSY 
corrections to pure QCD vertices. Actually, to have full control on the quantitative size 
of these effects we have explicitly checked that their numerical contribution is immaterial. 
To this end we have evaluated the full set of one-loop SUSY diagrams for the relevant 



tbH + vertex. The same set was considered in detail Ref. [10] in the case where all external 
particles are on-shell. In the present case, however, at least one of the quarks in that vertex 
is off shell (see Fig. 1). Therefore, we can use the same bunch of diagrams as in the on-shell 
case but we have to account for the off-shell external lines, which is a non-trivial task. We 
have studied this issue in detail by expanding the off-shell propagators 



1 



(p + q) 2 — M 2 q 2 - M 2 



1 _ p 2 + 2pq / p 2 + 2pq s _ 



2 -M 2 \q 2 -M 2 



(4.11) 



in the expression of the one- loop functions. The result is that the off-shell corrections 
are generally small for they are of order or below the non-leading one-loop effects that 
we have neglected. For the computation of the SUSY-corrected matrix elements we have 
proceeded in the following way: first, we have modified CompHEP's Feynman rules to 
allow for the most general ibH + vertex; then we have let CompHEP reckon the squared 
matrix elements and dump the result into REDUCE code. At this point, we have inserted 
expressions for the coefficients of the off-shell ibH + vertex that include the one-loop off-shell 
supersymmetric corrections to the vertex itself and to the off-shell fermion propagators and 
fermionic external lines 10 . Only half the renormalization of an internal fermion line has 
to be included, the other half being associated to the gqq vertex. Higher order terms not 



included in the effective Lagrangian (4.1) have been discarded. The subsequent numerical 
computation of the one- loop integrals has been done using the package LoopTools |{73|| . 
This procedure has allowed us to estimate the relative size of the off-shell effects in the 
signal cross-section, whose effects are analyzed in the next section. 

The upshot is that the approximation of neglecting qqg vertex and propagator cor- 
rections in the cross-section computation of process ( |2.1[ ), which may be called "improved 
Born" approximation (see next section), is really justified in the relevant region of param- 
eter space. 

4.2 Numerical analysis of the SUSY corrections 

We have seen above that the SUSY corrections are potentially large and can be of both 
signs. The leading correction to the cross-section of our process (|2.1| ) goes roughly as 

10 We shall not write down here the analytic expressions for the renormalized vertex and propagators. 
They can easily be derived by just generalizing previous on-shell calculations, such as those for t — ► bH + [hoi. 
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a — > a(l — 2AmJ at one loop, so that the sign of the SUSY-QCD part carries just the sign 
of —fx whereas the SUSY-EW one depends on the sign of the combined parameter — fxAf. 
Recall that the radiative B-meson decays (based on the b —* s 7 transition) prefer the sign 
fxA t < (see e.g Refs. |2(], |74|) if the charged Higgs boson is not exceedingly heavy. The 
reason is that if stops and charginos are of similar mass to that of the Higgs boson then 
they can compensate for the Higgs boson contribution, which by itself would overshoot 
the experimentally allowed range for BR(b — ► S7). In this case the fxAt < sector of 
the MSSM parameter space offers the Tevatron II a chance to observe a supersymmetric 
charged Higgs boson of around 200 GeV. For a much heavier Higgs boson, however, namely 
that accessible only to LHC, b — > s'y does not place any restriction on the sign of fx At. 11 

Obviously there is a preferred sign for our process: the optimal situation for the charged 
Higgs boson production mechanism ( |2.1| ) to be maximally sensitive to the SUSY effects is 
when fx < and A t > corresponding to both corrections, SUSY-QCD and SUSY-EW, 
giving a positive contribution to the cross-section. The next-to-interesting case is when 
/x < and A t < 0. Here the electroweak loops try to balance the SUSY-QCD ones, but the 
dominance of the latter still leaves a sizeable net outcome - some 60% of the previous case. 
Of course there is also the odd possibility that fx > and At > 0, but in this circumstance 
the SUSY quantum effects could still manifest themselves (if tan/3 is large) through an 
effective if -factor ( 2.1 2| ) sensibly smaller than the one predicted by QCD expectations: 
K MSSM < K^ CD . The worse possible situation would occur if tan [3 is small, say in the 
intermediate range 2 < tan (3 < 20, since no positive nor negative effect whatsoever could 
be detected; in fact, in this case not even the tree-level signal would be available! In the 
following we are going to concentrate our numerical analysis on the case fx < and At > 0, 
which is the one phenomenologically most appealing and therefore defines the scope of the 
present study. 

In Fig. H] we present the total cross-section as a function of tan/3 at the LHC (with 
M H + = 500 GeV) and the Tevatron (with M H + = 250 GeV) for K QCB = 1. We show the 
tree-level result and four additional cases that include SUSY corrections: they correspond 
to the four sets of parameters displayed in Table ||. The SUSY-corrected curves contain 
all the known corrections discussed in the previous section, including the resummation of 
the SUSY effects. As can be seen the numerical effect of the SUSY contributions can be 
dramatic; for the Set A there is a 100% positive enhancement of the cross-section. However 
the fact that we have resummed all the leading terms of the form (atan/3) n ensures that 
the prediction is robust under the inclusion of additional higher order effects. 

Let us comment briefly on the treatment of the strong structure constant as- We have 
used the SM value of as (Q) for the gluon-gluon and quark-gluon vertices. We have to do 
so, since the value of as is related to the determination of the PDFs and moreover the 
SUSY effects on these vertices are, as we have discussed before, non-leading in our case. 
However, in order to compute the SUSY corrections, we need to compute as(MsusY), and 
it would be a very bad approximation to use the SM value. So we use the MSSM value 



11 The precise measurement of the anomalous magnetic moment of the muon j75| could favor /i > [ [T6[ , 
but the impact on our case is marginal because the (g — 2) M value in the MSSM critically depends on the 
values of some slepton masses whose influence on our calculation is completely negligible. 
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Table 6: Sets of SUSY parameters used in the computation of the SUSY-corrections to the charged 
Higgs boson associated production (all masses and trilinear couplings in GeV). 
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Figure 9: Total cross-section for the process pp(pp) — ► H + ib + X as a function of tan/3 in the 
relevant range of this parameter. The tree-level result and the SUSY-corrected result are shown 
(with K^ CT> = 1) for the four sets of parameters from Table ||, for a) the Tevatron with M H + = 
250 GeV, and b) the LHC with M H + = 500 GeV. 



a^ SSM (Q) to compute the corrections to the H + ib vertex. 

Of course, one would like to know which is the relative importance of each effect in 
the final corrected cross-section for the process fl2.1|). Let us compare the following set of 
corrections: 

• the full set of SUSY corrections, including resummation and the off-shell two- and 
three-point functions (labeled Full); 

• the corrections including resummation, but using on-shell values for the irreducible 
three-point functions {On- Shell); 

• the one-loop (without resummation) corrections, using also on-shell values for the 
irreducible three-point functions (OS-nr); 

• the tree-level result with the pole quark masses replaced by the effective Yukawa 



couplings eq. (4.1) (Improved- Born). 



In Fig. |l0|we show the total cross-section of process (|2.l| ) using the four approximations. 
We use the intermediate Set B of Table which represents a moderate case. The following 
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Figure 10: As in Fig. 9, but now showing the different approximations described in the text for 
the intermediate Set B in Table ||, a) for the Tevatron with M H + = 250 GeV, and b) for the LHC 
with M H + = 500 GeV. 



hierarchy of quantum effects is observed: 

a(Full) > ^(Improved - Born) > a(On - shell) > cr(OS - nr) . (4.12) 

The effect of the resummation is clearly visible. The difference between the resummed 
and non-resummed results grows with the magnitude of the corrections; in sets A and D 
this difference is much larger than that of sets B and C. On the other hand the difference 
between the Full and the On-Shell results is nearly indistinguishable. So the Off-Shellness 
effect is generally small in our case and it may be neglected in a first approximation, 
although it should be taken into account in a detailed study. 

In the following we restrict ourselves to plot the effective SUSY K-factor (K slJ ) 
involved in the definition ( 2.1 2| ) . This factor embodies the main results of our calculation 



as it gauges directly the potential impact of the genuine SUSY corrections. It is defined as 
follows: 

_.SUSY-corrected 

*- SUSY = jT „.-,„, ■ <«»> 



for the different contributions mentioned above. 12 In Fig. 11 we present the evolution of the 



i^ SUSY -factor ( |4.13 ) as a function of various MSSM parameters for the case of the Tevatron 



and for a charged Higgs boson mass of 250 GeV. In Fig. ^ we show the same effective 
if SUSY -factor for the LHC and for a charged Higgs boson mass of 500 GeV. Both cases 
look pretty similar, the reason being that the leading corrections (eqs. ([Q-4.5), Fig. |8|) 



are insensitive to the charged Higgs boson mass and the energy of the process. For the 



12 The pole quark masses are used in the ratio (4.13). Notice that although the total cross-section changes 
significantly by the choice of running or on-shell masses, the change in K SVSY j s rather mild, since it affects 



both the numerator and denominator of (4.13) 
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Figure 11: if SUSY -factor - eq. ( [4.13D - for the parameter Set B -Table |- and M H + = 250 GeV at 
the Tevatron, as a function of a) the gluino mass; b) the /i parameter; c) the lightest stop mass and 
d) the top-squark soft-SUSY-breaking trilinear parameter A t . Shown are the four approximations 
explained in the text. 



same reason the main features of the corrections follow closely the pattern already observed 
in the partial decay widths T(t — ► H + b) [JTo|] and T(H + — > tb) In these figures we 

have concentrated on the moderate Set B, although larger squark masses will not decrease 
significantly the corrections, unless there exists a large hierarchy between the masses of the 
squarks, the gluino and the [i parameter. 

For the parameter Set B, the corrections are always positive. They decouple very slowly 
with the gluino mass, as seen in Figs. [Tl] a and [l^a, so a very heavy gluino does not prevent 
to have large corrections. In Figs, [ll] and 12 the role of the resummation is very clear; failing 
to include it would lead to completely underestimated cross-sections. The shaded regions 



in Figs, lib and |l2Jb correspond to a chargino mass below the LEPII mass exclusion limit. 
The parameter is also a key parameter in the process under study, since the SUSY-QCD 
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Figure 12: JT SUSY -factor - eq. fl4.13p - for the parameter Set B -Table ||- and M H + = 500 GeV 
at the LHC, as a function of a) the gluino mass; b) the /i parameter; c) the lightest stop mass; and 
d) the top-squark soft-SUSY-breaking trilinear parameter A t . Shown are the four approximations 
explained in the text. 



corrections grow linearly with it - Cf. Eq. (f4.3|). In the \i < scenario the SUSY-QCD 
corrections to the cross-section are positive: roughly, 8a jo ~ — 2 (Arrib)susY-QCD > 0. 
Notice, however, that in Figs, lib and |i~2|b, At > 0, and so the SUSY-EW corrections add 
constructively to the SUSY-QCD ones bo jo ~ — 2 (Amb)suSY-EW > 0. In the /i < 
and At < scenario, instead, the SUSY-EW corrections would partially compensate the 
positive SUSY-QCD loops, giving a if SUSY factor slightly smaller. Finally, there is the 
case n > and A t > 0, where the two SUSY effects on the cross-section would be negative: 
if SUSY < 1. As we have already mentioned, this could lead to if MSSM ~ 1 in eq . ( 2 . 1 2| ) 
and then the signature of the underlying SUSY would be that the QCD corrections are 
"missing" or even negative! While this is possible for an Htt final state p7| 1, for example 
due to Coulomb exchange of gluons between quasistatic top quarks, this is not allowed in 
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(a) 



(b) 



Figure 13: if SUSY -factor - eq. (4.13) - for the parameter Set B -Table g- for the different 
production channels as a function of tan/3: a) For the Tevatron, and M H + = 250 GcV and b) the 
LHC and M H + = 500 GeV. 



the presence of 6-quark in final states at high energies. 

We come now to the effects from the squark sector. In Figs. 11c and 12c we show 
the evolution with the corrections with the lightest stop-quark mass. The plots end at 
TO£ = 500 GeV because for larger stop-quark masses there is no physical solution for the 
mixing angle (given = 500 GeV) in this particular set of inputs. We see that, although 
the corrections have a dependence on the squark masses, this parameter is not the most 
critical one. In fact, we emphasize that larger squark masses (as e.g. mg — ~ 1 TeV) 
do not necessarily lead to smaller corrections - see Fig. ^, parameter Set A. 

We finally focus on the influence of the soft-SUSY-breaking trilinear stop-quark cou- 
pling At (Cf. Figs, lid and p"2]d), which is another key parameter for the corrections under 
study. We see that, again, the preferred range (// < 0, A t > 0) gives the largest positive 
corrections, thus enhancing the total cross-section dramatically. 

The final point that we wish to address is the study of the relative contributions from 
the different channels involved in the production cross-section -eqs. ( |2.2| )-( |2"~4|) . Since the 
main effect comes from the Araj term, we expect that the SUSY-corrections will be very 
similar in all the channels. In Fig. [D| we display the corrections from the different channels 
for the Tevatron and the LHC. We see, indeed, that the corrections to the various channels 
are very similar. This means that to obtain the discovery potential (or exclusion region) 
for the different accelerators in the presence of the SUSY-corrections, we may directly scale 
up the previous results from Fig. [7| by a factor K SUSY without having to remake the full 
kinematical analysis. 



4.3 Consequences for the charged Higgs boson search 

In the previous figures we have exemplified a wide variety of MSSM effects on the 
cross-section of process ( [2.1D based on a typical set of MSSM inputs (sparticle masses and 
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soft-SUSY-breaking parameters). The new contributions entail cross-section enhancements 
at the level of 50 — 70% beyond standard QCD expectations. In other words, the SUSY 
loops have the capability to modify the largest expected QCD A'-factor from AmS° = 1.6 
up to an effective MSSM value (pl^ ) reaching K MSSM = 2 A- 2.7. Surely enough, if such a 



dramatic effect is really there, it could not be experimentally missed. In turn the discovery 
and exclusion limits on the charged Higgs boson inferred from the previous analysis will 
get significantly modified as compared to the tree- level case (Cf. Fig. ^), as we shall see 
next. 



As already warned in Section p-lL one should also keep in mind the SUSY corrections 
to the branching ratio BR(H + — > tb). They are certainly much tamed than those poten- 
tially affecting the partial decay width T(H + — > tb) j77|], but even so they can have some 
impact on the charged Higgs boson mass discovery /exclusion limits. Recall that the QCD 
corrections decrease significantly the partial decay width T(H + — > tb) f77j| . The genuine 
SUSY corrections, on the other hand, are governed by the same kind of contributions dis- 



cussed in this paper, and can therefore have both signs [12, 35|. As the sign of the SUSY 
quantum effects on the branching ratio is the same as for the production cross-section, they 
just mutually enhance or suppress the signal. Here we use the improved expression for this 



partial decay width given in Ref. [35|. The only other decay channel that we consider is 



H+ — ► t + u t , but we neglect the SUSY corrections to its partial decay width. Notice that, 
since BR(H + tb) is much larger than BR(H + — > t + v t ), the quantum corrections to the 
branching ratio BR(H + — > tb) will be small in general. We remark that in the light charged 
Higgs boson mass scenario relevant to our tibb signature (viz. mt + m^ < M H ± < 300 GeV), 
the partial decay channel H + — > t + v t has a sizeable branching ratio which raises very fast 
with decreasing M H ± . For instance, for M H ± ~ (230 — 240) GeV and a moderately high 
tan/3 = 30 one gets BR{H + — > t + v t ) > 20% (without SUSY corrections). In some cases 
(specifically when /u > 0, At < 0, like in Sets C and D in Table |6|) this branching ratio can 
be further enhanced by SUSY effects due to the corresponding decrease of F(H + — > t + b) 13 . 
Thus for charged Higgs masses in the heavy range (Mjj± > 300 GeV), where the uncor- 
rected branching ratio of the r-mode would be below 10%, the latter can actually remain 
at the level of 15% (even for M H ± > 500 GeV) thanks to the SUSY effects; this scenario 
would be realized e.g. for {m~ g ,m ix , m~ h , n, A t ,A b ) ~ (500, 150, 300, 200, -200, 200) GeV at 
tan/? = 30 as before. Hence, it is clear that there are situations for M H ± > rat — "ifc where 
the r-mode could be successfully used to explore the light charged Higgs mass region, as 
it was already the case in Run I for M H ± < mt — rrih through the study of the combined 
decay t — > H + + b — > t + u t + b [|lj], |1| 16, 17]. Nevertheless, the r-mode can alternatively 



be highly suppressed by the SUSY effects themselves - if they turn out to greatly enhance 
T(H + — > t + b). Indeed, for (m^ , , ) ~ (500, 150, 500) GeV and assuming (fi,At,Aj,) 
as in Sets A and B, the r-decay of the charged Higgs narrows down to a 5% or less. There- 
fore, the SUSY suppression of the r-mode turns out to occur in the regions of the MSSM 
parameter space where the production processes pp(pp) — > H + tb + X are maximized. We 
conclude that in the parameter space relevant for pp(pp) — > H + tb + X the tibb signature 



13 See J.A. Coarasa et al. in fij) for the SUSY analysis of BR(H + -» r + u T ). 
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is in general the most suited one, and indeed it is the only one used in the present study. 
We warn the reader that at the opposite end of the spectrum, where M H + > 1 TeV, the 
branching ratio BR(H + — * tb) would significantly decrease if M H + > + , due to 
the opening of the squark decay channels H + — ► b$j - which quickly become dominant 
at high tan/3. Given the parameters in Table [6|, these squark channels are relevant only 
for Mjj+ > ITeV (resp. 2 TeV) for the parameter Sets B and C (resp. A and D). Since 
these squark channels have not been considered in the present study, the conclusions on the 
charged Higgs boson discovery /exclusion limits are strictly valid only for M H + < +mf 1 . 
While we could easily include the new H + decay modes, the typical results from our study 
would not change at all. In fact, we have emphasized that the leading SUSY corrections 
relevant for our process do not necessarily decrease with the SUSY mass scale, so that 
our numerical results can be extended to > 1 TeV since we can always find a mod- 

ified particle spectrum satisfying M H + < + and producing the same radiative 
corrections as the corresponding parameter set in Table ^. 

In Fig. 14 we show the equivalent of Fig. [7], but using the full SUSY-corrected cross- 
section, and the corrected branching ratios, for tan/? = 50. In the Tevatron case, whose 
tree-level cross-section borders the discovery limits from below, we have applied a QCD 
correction factor Jf Q CD = 1.5 to account for the (yet unknown) standard QCD corrections 
to the production cross-section. Without large QCD corrections, the Tevatron cannot dis- 
cover a heavy charged Higgs boson M H + > 200 GeV (not even at a 3a level of significance) 
unless the SUSY corrections are as high as e.g. those corresponding to Set A in Table [| In 
this case a charged Higgs boson in the range M H + < 255 GeV can still be discovered at 3a 
at the Tevatron. This is certainly a welcome possibility as it would throw light on the SUSY 
nature of that charged Higgs boson, in case of being detected, and in addition it would 
strongly hint at regions of the MSSM parameter space very much akin to those defined by 
Set A. On the other hand, for intermediate SUSY sets (like Set B), the concurrence of the 
QCD corrections is indispensable to get a similar charged Higgs boson signal. At the LHC, 
in contrast, large discovery/exclusion regions exist already using the tree-level approxima- 
tion 

(i^QCD = 1); 

so we are more conservative in this case and we can afford not to apply 
a if Q CD correction factor. The band in Fig. [l4| interpolates the parameter space between 
the parameter sets A and D -Table |^-, which can be regarded as plausible parameter sets 
defining the typical range of maximum and minimum cross-sections expected in the MSSM. 
The interpretation of this figure is as follows: in the context of the MSSM, and within the 



representative region of parameter space discussed in Section 4.1, the Tevatron should be 
able to discover a relatively heavy charged Higgs boson {Mjj+ < 250 GeV) at the canonical 
5(7 significance only in an scenario of very large positive SUSY corrections. Nonetheless, 
the Tevatron can still find some reasonable evidence (at 3a, i.e. 99.7% C.L.) of a charged 
Higgs boson with a mass 200 — 250 GeV if only the condition fj, < applies irrespective of 
the sign of At- The LHC, on the other hand, does guarantee a 5a discovery in the range 
M H + < 800 GeV for 100 fb _1 of integrated luminosity; and a 3a evidence of a charged 
Higgs boson with a mass up to 1 TeV for any SUSY spectrum. With only 30 fbT 1 the LHC 
can not discover a superheavy charged Higgs boson M H + > 1.1 TeV, but the discovery is 
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Figure 14: Total effective cross-section for the process (2.1) as a function of the charged Higgs 
boson mass a) at the Tevatron with a QCD correction factor i^Q CD = 1.5 and b) the LHC (with 
^QCD = i). In both cases tan/3 = 50. Shown is the total cross-section interpolating between 
parameters sets A and D (shaded region) and for the parameter Set B. Also shown are the cross- 
sections necessary for discovery and exclusion from Fig. [?]. 



possible for M H + < 1.1 TeV and fi < 0. Once we will have data from the LHC this picture 
could change significantly. If SUSY particles are found, and we know (roughly) their mass 
scale, the band in Fig. 14 will get significantly shrunk, providing clean predictions. If, 
on the contrary, a charged Higgs boson is found, but the scale of the SUSY particles is 
still not known with sufficient accuracy, one can use Fig. 14 to extract information on the 
SUSY particles properties. Clearly the potential impact of the SUSY corrections is not at 
all negligible and should not be ignored in any serious analysis of this subject. 

By iterating this procedure we find the regions of the tan/3 — Mjj+ plane in which each 
experiment can find (or exclude) the existence of the charged Higgs boson. We present 
the result of this analysis in Figs. p~5| a and b. These figures make transparent what was 
stated above: the presence of the SUSY corrections alters significantly the Higgs boson 
discovery potential of the hadron colliders. From the figures it is also patent the Tevatron 
could have a non-negligible chance to find a charged Higgs boson in the intermediate range 
M H ± < 280 GeV. 



We have also shown in Fig. 15 the value of the pseudoscalar Higgs boson mass M^o 
corresponding to each value of M H + . We have correlated M^o with M H + using the one-loop 
corrections to the Higgs bosons masses as provided by the program FeynHiggsFast ||78|| . 
Note that in this case Mjj+ changes with tan/?, but is also different for each set of SUSY 
parameters. However, the change of M^o among the different sets is very small (below 
5%), and the one-to-one correspondence M H + «-> M A o is a very good approximation. 

Since the performed analysis of the SUSY corrections is meaningful only for large tan/? 
(tan/3 > 20), the shifts in the discovery regions in Fig. 15 b are not sufficiently reliable for 
tan/? < 20. Unfortunately this means we can not compare with the results of Ref. 0], 
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Figure 15: Discovery /Exclusion regions for the Tevatron and the LHC, in the tan/3 — M H + plane. 
Also shown is the value of M^a corresponding to each M H + according to the MSSM. 



given for M^o < 500 GeV (which correspond to tan/3 values < 20 in Fig. 15b). Sticking 
to the large tan/3 scenario, we see that for tan/3 = 60 a charged Higgs boson discovery is 
guaranteed for M H + < 1 TeV. For /i < a charged Higgs boson could be discovered up to 
a mass of 1 TeV for any tan/3 > 30. 

We should also recall here that our results show a certain degree of uncertainty due to 
the dependence of the cross-sections on the particular PDF parametrization and the vari- 
ation with the renormalization and factorization scales, which we have already illustrated 
in section ||. These issues are particularly relevant for the Tevatron, where a charged Higgs 
boson of order 250 GeV could be discovered only if the QCD and SUSY effects add up 
with the same positive sign. In the worst possible scenario, the K^ CY) = 1.5 factor that is 
needed in combination with the typical SUSY effects from the parameter Set B in Table || 
could be effectively reduced. In this adverse circumstance the charged Higgs discovery at 
the Tevatron would only be possible if the SUSY effects would be still larger, namely at 
least as large as those from the parameter Set A in that table. In any case the lesson is 
once again the same: without the collaboration of SUSY (or some other form of physics 
beyond the SM), the charged Higgs discovery is impossible at the Tevatron. On the other 
hand, the aforementioned uncertainties are much more tamed for the LHC, and our results 
shown in Fig. [l4|b and 15b should not be too much affected by them. 



5. Summary and conclusions 

We have presented a thorough study of the process pp(pp) — > H + ib + X in hadron colliders 
in order to assess the possibility to see a SUSY charged Higgs boson at the Tevatron II and 
at the LHC. A necessary starting condition for this is that tan/? is sufficiently large (viz. 
tan/3 > 20) so as to insure that the cross-section of the process pp(pp) — > H + tb + X is 
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sizeable. Beyond this tree-level requirement, our study of the quantum corrections to that 
process within the MSSM has shown that they are dominated by exceptionally important 
effects that can be absorbed into an effective tbH + vertex, and therefore in practice they 
can be treated at the level of an "improved Born approximation". After the inclusion 
of the MSSM effects we have found that the "effective K-factor" ( p. 12 ) for the process 
Pp(pp) ~~ * H + ib + X could be either smaller than expected, perhaps barely if MSSM < l } or 
on the contrary substantially higher, sayif MSSM = 2-3. In both cases it should be possible 
to exploit this process to unravel the nature of the new physics in which the purportedly 
discovered charged Higgs boson is integrated. This could be performed by detecting the 
existence of large radiative corrections beyond the SM, namely sizeable quantum effects 
that remain after subtracting the full set of conventional QCD corrections at the NLO. This 
will typically be the case if the squark masses are heavy at an intermediate level of order of 
a few hundred GeV. The gluino, notwithstanding, can be as heavy as 1 TeV. Remarkably 
enough, in certain circumstances explained in the text all the sparticle masses, and not 
only the gluino, could be of order of 1 TeV, and still the MSSM quantum effects would 
remain large (~ 50%). This feature shows that sound hints of SUSY can be derived from 
the sole study of Higgs boson production even if the whole sparticle spectrum is barely 
accessible to the future colliders. 

If tan/3 > mt/mb > 35, then the upgraded 2 TeV-Tevatron can exclude charged Higgs 
boson masses below 300 GeV at the 95% C.L. On the other hand, it is not likely that 
the Tevatron will be able to produce any signal at the 5<7 significance level to "officially" 
discover a SUSY charged Higgs boson. The only chance would be in the narrow window 
M H + = (180 — 250) GeV and assuming the most favorable region in the MSSM parameter 
space (typically, Set A in Table ||). In any case, the Tevatron can provide reasonable 
evidence (at the 3 a level) of a M H + = (200 — 300) GeV Higgs boson, and also it can do a 
good complementary job at high tan (3 using the neutral Higgs bosons channels pp(pp) — > 
hbb + X (h = h , H°, A ) already studied in the literature B |3L ffl, where the absence of a 
top quark in the final state increases the signal in some of them provided the Higgs boson 
is not too heavy. In this respect we have seen that in the light charged Higgs boson mass 
scenario, the partial decay channel H + — ► r + v T has a large branching ratio, and in some 
favorable cases it could be successfully used at the Tevatron II to explore this region, as it 
was the case for the Run I [14, 15, 16, |l7|]. Fortunately, the situation with the LHC is much 
more rewarding under similar circumstances. In fact, the LHC will be able to discover a 
charged Higgs boson up to 800 GeV or else to exclude it up to 1 TeV at least. Remarkably, 
these limits can be significantly improved up to 1.1 TeV and 1.5 TeV respectively in the 
presence of SUSY corrections. 

At the end of the day the strategy should be clear and can be put in a nutshell. If 
the signal pp(pp) — * H + tb + X is found it will automatically imply new physics. Moreover, 
under the assumption that the discovered charged Higgs boson is a member of a Type II 
two-Higgs-doublet model it will strongly suggest a high value of tan /?, certainly not below 
20, and most likely above 30. Then an accurate measurement of the tbH + vertex should 
be possible at a level better than 10%, together with its correlation with the neutral Higgs 
bosons vertices bbh {h = h°,H°,A°) obtained both at the LHC and at the Tevatron from 
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the measurement of the companion processes pp(pp) — > hbb + X. After subtraction of the 
conventional QCD corrections the genuine non-standard quantum effects could possibly be 
disentangled after accurate comparison of all these signals [31, |3(J. As we have shown, the 
"remnant" of this subtraction could be quite large in the MSSM case, viz. of the order of 
the full QCD effects themselves! Then a pattern of sparticle masses (or at least of favoured 
regions in the MSSM parameter space) responsible for these virtual effects could emerge 
(see e.g. figures p|-|l2|), and of course it should be cross-checked with the results obtained 
from other independent experiments, e.g. at the LC j?^]. Altogether the bottom line of 
this strategy could definitely point towards the supersymmetric nature of the found Higgs 
bosons until the sparticles themselves will eventually be produced and correctly identified. 
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